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Abstract We perform a lattice study of the topological sus- 
ceptibility and instanton size distribution of the SU(2) gauge 
theory with two adjoint Dirac fermions (also known as Mini- 
mal Walking Technicolor), which is known to be in the con- 
formal window. In the theory deformed with a small mass 
term, by drawing a comparison with the pure gauge theory, 
we find that topological observables are decoupled from the 
fermion dynamics. This provides further evidence for the in- 
frared conformality of the theory. A study of the instanton 
size distribution shows that this quantity can be used to de- 
tect the onset of finite size effects. 
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1 Introduction 

After the recent experimental breakthroughs that allowed 
the announcement of the discovery of a new bosonic par- 
ticle with mass at the electroweak scale IDE], the focus of 
experiments is moving towards the determination of the na- 
ture of this boson, with the goal of finding hints for new 
physics. Hence, it is of utmost importance for theorists to in- 
vestigate possible scenarios of electroweak symmetry break- 
ing beyond the standard model. An appealing possibility is 
the existence of a new strong interaction whose chiral con- 
densate breaks electroweak symmetry. Based on this idea, 
over the years a physically consistent framework denoted 
as Walking Techinicolor (WT) was developed 1 3 1 4 , 5 ,|6]|7][8l 
HnO||n]|ll| (reviews on the topic include (BEE)). In 
order to be consistent with the stringent electroweak preci- 
sion measurements, a theory realising the WT scenario must 
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have two crucial properties: dynamics characterised by near- 
conformality, and a chiral condensate with anomalous di- 
mension of order one. In most semi-analytical treatments of 
WT, both properties are taken as reasonable assumptions. In 
order to prove that there exist realistic gauge theories char- 
acterised by near-conformal behaviour and large anomalous 
dimension, calculations from first principles are needed. Gauge- 
string duality techniques can be used to engineer theories 
with the required properties in the limit in which the number 
of gauge bosons goes to infinity (see e.g. lll6lll7l[T8l[T9ll20l 
|2D and references therein). For a gauge theory based on a 
finite Lie group, lattice computations provide the best quan- 
titative tool. 

In recent years, a large body of activity has been devoted 
by the lattice community to computations aimed at charac- 
terising the phases of gauge theories, with the goal of un- 
derstanding signatures of the conformal and near-conformal 
regime. The current understanding is summarised in ll22l . 
Among the theories investigated to date, SU(2) with two 
Dirac fermions in the adjoint representation (also known 
as Minimal Walking Technicolor [23,24,25], for which ex- 
perimental constraints have been recently reported in |26|) 
plays an important role: it is the only theory studied so far 
to show unambiguously the expected behaviour of an in- 
frared conformal theory, in independent numerical results 
obtained with different techniques and targeting different 
observables. Although for building a technicolor model one 
needs a near-conformal gauge theory, first principle stud- 
ies of infrared conformal theories are interesting both from 
a theoretical point of view and from a phenomenological 
perspective GTll . After the first exploratory calculations of 
physical quantities in this theory [28 29 30 31 1, which were 
fundamental in ascertaining its conformal behaviour, lattice 
studies have focussed on the running of the coupling near the 
fixed point (32,33,34, 35 1, on the calculation of the anoma- 
lous dimension of the chiral condensate [36, 37, 38, 3911401 
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I4T1 and more recently on controlling finite size effects 11421 
I43ll44l . The main goals of this community effort have been 
to establish a reliable technique to extract the anomalous di- 
mension and to identify a good set of observables that could 
provide an unambiguous signature of conformal behaviour. 
These investigations are expected to be relevant also for the 
near-conformal case, in which the theory will mimic the 
properties of a conformal theory at some intermediate dis- 
tances before showing its true confining and chiral symme- 
try breaking nature at asymptotically large distances. 

To date, the attention of lattice studies has focussed on 
spectral quantities. In this paper, we shall report on a first 
investigation of topological quantities^] and in particular on 
the behaviour of the topological charge, of the instanton size 
distribution and of the topological susceptibility. These quan- 
tities play an important role in confining and chiral symme- 
try breaking theories like QCD (see e.g. flTl for a recent 
work), but their behaviour in (near-)conformal gauge theo- 
ries has not been explored before. We shall study whether 
those observables can be used to identify the infinite volume 
regime and the scaling regime in the mass-deformed the- 
ory. Taking inspiration from the perturbative argument given 
in B8ll . in |37| it has been shown from first principles that 
SU(2) gauge theory with two Dirac adjoint flavours behaves 
like a confining theory with heavy adjoint quarks for any fi- 
nite value of the fermion mass, which plays the role of a de- 
forming parameter in a would-be infrared conformal theory. 
Hence, in the infinite volume limit, we expect topological 
observables to be compatible with those of Yang-Mills, with 
any small deviation originating from the fermionic matter; 
large deviations are due to finite size effects. In addition, 
in the large volume and small mass regime, the topological 
susceptibility % will show the scaling behaviour 



X °= m l +f 



(1) 



with y* the anomalous dimension of the chiral condensate 
and m the fermion mass. At finite but large volume L 4 , this 
relationship becomes 

X =L- 4 F(L l+ ?*m), (2) 

with L the linear size of the lattice and F a function of the 
scaling variable x — L m. 

The purpose of this investigation is twofold. On the one 
hand, we aim to verify that for large volumes topological 
variables are broadly compatible with Yang-Mills values and 
at characterising qualitatively their behaviour on smaller vol- 
umes. On the other hand, we want to investigate whether the 
scaling relations ([T} and Q can be used for an accurate mea- 
surement of the anomalous dimension. The rest of the paper 
is organised as follows. In Sect. [2] we review the formula- 
tion of the theory of interest on the lattice and we define the 
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observables related to topological quantities that are studied 
in this work. Our numerical results are reported in Sect. [3] 
Finally, Sect.[4]summarises the main findings of our study. 



2 Methodology 

In this work, we study the theory on an four-dimensional 
Euclidean spacetime lattice. The lattice discretisation used 
here is described in [29 1 . We recall the main points below, 
referring to |29| for further details. 

On the lattice, the dynamics of the gauge variables is de- 
scribed by group elements in the fundamental representation 
living on the links of the lattice. If (i) is the field living on 
the link originating from the lattice point i in the direction 
fi, the plaquette operator t/^ v (0 ls defined as 

EWO - Up(i)Uv(i + fL)U*(i+ v)ul(i) . (3) 

In SU(N), the Wilson action for the gauge degrees of free- 
dom of the theory is then 

1 
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with j3 = 1 /g 2 and g the gauge coupling. In the Wilson for- 
mulation, which is the discretisation used in our work, the 
dynamics of the fermions is described by the Dirac operator 
M a p(ij), defined as 

M a p(ij) = (m + 4r)Sjj8 a p - ~ 



hJ+H 



(5) 



where i and j are lattice sites and a and j5 Dirac indices. 
The superscript R indicates that the link variable have to 
be taken in the representation of the fermions and m is the 
bare fermion mass. Since this fermion discretisation breaks 
chiral symmetry, the chiral point has to be determined non- 
perturbatively in the numerical simulation. A possible def- 
inition of this point uses the fermion mass mpcAC defined 
through the partially conserved axial current: the chiral limit 
is attained at the point for which mpcAC is zero. 

With those definitions, the lattice discretised version of 
the SU(2) gauge theory with two adjoint fermions we are 
studying is described by the path integral 



Z = J (W M (0) (detM(f/ M )) w / e - s 



(6) 



where the link variables are SU(2) matrices, M is in the ad- 
joint representation and the number of flavours Nf is 2. A 
set of configurations that dominate the path integral can be 
generated with Monte Carlo methods, and vacuum expec- 
tation values of observables are computed as ensemble av- 
erages of the corresponding operators over those configu- 
rations. These observables have a relative statistical fluctu- 
ation that scales as 1 / t/Ns, where N s is the number of the 
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generated configurations. Thus the numerical error can be 
kept systematically under control by adjusting the size of 
the sample. 

The observables we have studied are described in [49 ], to 
which we refer for further details about the observables and 
the model assumptions made to derive the relationships we 
use in our work. Here we only summarise the main points. 
In the continuum, the topological charge density can be ex- 
pressed as 



GW 



l 



Tr{F MV (x)F pt7 (x)} 



(7) 



32n 2 C " vp 

The total topological charge of a configuration can then be 
obtained as the integral of this quantity over the spacetime. 

The equivalent lattice topological charge density is then 
given by: 



Gl(0 
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e M vpffTr{t/ |lv (/)J/ p(T (i)} • 



(8) 



For a smooth gauge field, this would have fluctuations 
of order a 2 as the lattice spacing a is sent to zero; how- 
ever, realistic fields have ultraviolet fluctuations which will 
completely dominate over the physics of interest in the con- 
tinuum limit. To mitigate this, a cooling process is intro- 
duced |50|. The cooling process operates by minimising the 
local action for each lattice site in turn. Successive cool- 
ing sweeps will "smooth out" the fluctuations such that the 
physics may be observed. A side effect however is that in- 
stantons may be shrunk or may annihilate in instanton-anti- 
instanton pairs, thus excessive cooling is to be avoided. The 
number of cooling sweeps performed in a calculation results 
from the compromise between the need to smooth out the 
configurations and the necessity of not losing physical in- 
stantons by annihilation. The tuning of this parameter is not 
critical, since there is a wide plateau where the physics can 
be observed before cooling artefacts set in. 

Once the topological charge density is known, various 
observables may be calculated. The topological susceptibil- 
ity is defined as 
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(9) 



where V is the lattice volume and Qj is the total topological 
charge, defined as 



Gr = £fiL(0 



(10) 



with i running over all lattice points. Since in the continuum 
the total topological charge is an integer, on the lattice Qj 
is often rounded to the nearest integer. Alternative discreti- 
sations can be used. Lattice studies (e.g. [51 52|) show that 
the particular definition of Qj does not affect the results if 
the calculation is done close to the continuum limit, as one 
would have expected. 
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Fig. 1 Variation of the quantity % 3 / \f® w ith o. 

The size of a given instanton may be calculated from the 
local maxima of the absolute value of the topological charge 
density given in Eq. d8) from the relation 



<7peak 



n 2 p 4 



(11) 



which has been used to determine the instanton size distri- 
bution and the average instanton size. 



3 Results 

Preexisting two-flavour configurations discussed in |[39l (to 
which we refer for measurements of spectral quantities and 
for a careful discussion of finite volume artefacts) were used, 
at J3 = 2.25, with V = 16 x 8 3 , m = -0.90, -0.75, -0.50, 
-0.25, 0.00, 0.25, 0.50 and V = 24 x 12 3 , m = -1 .00, -0.95. 
20 cooling sweeps were used for all configurations. For com- 
parison with the Yang-Mills case, pure gauge configurations 
were generated at J3 = 2.4, V = 12 4 and j8 = 2.55, V = 20 4 ; 
each update consisted of 1 heat bath and 4 over-relaxation 
steps, and measurements were taken every 10th and 50th 
configuration on the smaller and larger lattice respectively. 
At each /3, the size of the lattice has been chosen so that 
finite size effects are negligible. 

We start with the topological susceptibility. For this quan- 
tity, we have used the integer definition of Qj, after having 
verified that alternative definitions give compatible results. 
Our findings (figure [TJ show consistency between the pure 
gauge and the 2-flavour theory, as one would have expected. 
As a check, we have verified that our pure gauge results are 
consistent with those presented in [51 1. Although no sign 
of finite size effects can be seen in the topological suscepti- 
bility, the instanton size distribution tells another story. We 
start by comparing this observable in the pure gauge case 
and in the two flavour theory on the larger lattice. The in- 
stanton size distribution is consistent in these cases (figure 
[2]), which in turn are consistent with the distribution found 
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Fig. 2 Instanton size distribution for the two-flavour configuration sets 
and for the pure gauge configuration sets at the values of the lattice 
parameters shown. 



Fig. 4 Comparison between the instanton size distributions for the dy- 
namical case at the smaller and larger lattice volume, for the sets of 
bare parameters shown. 
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Fig. 3 Instanton size distribution in the dynamical case on the smaller 
lattice at the parameters shown. The appearance of a peak at larger 
instanton size signals the onset of finite size artefacts. 
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Fig. 5 Variation of the average instanton size with a. 



in similar studies 1 53,54,55]. However, on the smaller lat- 
tice the Nf = 2 theory shows a long tail that develops onto a 
second peak as the region affected by small size artefacts is 
entered (figures[3]and[4]). As a consequence, when finite size 
artefacts are present the average instanton size looks larger 
(figure |5). 

Finally, we performed a fit of the topological suscepti- 
bility data according to Eqs. ([T]l and pi, but we did not see 
the expected scaling. The most plausible explanation is that 
larger lattices are required in order for the scaling behaviour 
to manifest. Generating configurations in the scaling region 
of gluonic observables requires a big computational effort, 
which is outside the scope of this work. Currently, a calcula- 
tion in this direction is being performed [44 1 . These configu- 
rations will allow us to investigate the scaling of the topolog- 
ical susceptibility as the mass goes to zero in the appropriate 
range of volumes and then to answer the question of whether 
the topological suceptibility is an efficient observable to ex- 
tract the anomalous dimension. 



4 Conclusions 

In this work, we performed the first numerical investigation 
of the behaviour of instanton-related quantities in a mass- 
deformed infrared conformal gauge theory. For the theory 
we have studied, SU(2) gauge theory with two dynamical 
adjoint Dirac fermions, we found that the behaviour of the 
topological susceptibility supports the scenario in which the 
infrared dynamics is dominated by gluonic quantities, with 
the fermions being more massive 1 48 . 37 1 . Moreover, the in- 
stanton size distribution provides a good indication of finite 
volume effects. We plan to extend this study to larger lat- 
tices, in order to check whether the topological susceptibil- 
ity obeys the expected scaling with the mass and whether 
we can extract from it the anomalous dimension of the con- 
densate. Finally, it would be interesting to apply those tech- 
niques to other theories that could provide a realisation of 
the walking scenario, like e.g. SU(3) gauge theory with two 
sextet fermions I56ll57l. 
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